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The effect induced by an environment on a composite quantum system is studied. The model 
considers the composite system as comprised by a subsystem A coupled to a subsystem B which 
is also coupled to an external environment. We study all possible four combinations of subsystems 
A and B made up with a harmonic oscillator and an upside down oscillator. We analyzed the 
decoherence suffered by subsystem A due to an effective environment composed by subsystem B 
and the external reservoir. In all the cases we found that subsystem A decoheres even though it 
interacts with the environment only through its sole coupling to B. However, the effectiveness of the 
diffusion depends on the unstable nature of subsystem A and B. Therefore, the role of this degree 
of freedom in the effective environment is analyzed in detail. 

PACS numbers: 03.65.Yz,03.65.-w,05.40.-a 



I. INTRODUCTION 

Decoherence is the process by which most pure states evolve into mixtures due to the interaction with 
an environment 0. The very notion of a quantum open system implies the appearance of dissipation and 
decoherence as an ubiquitous phenomena and plays important roles in different branches of physics |2| (from 
quantum field theory many body and molecular physics to theory of quantum information), biology 
and chemistry. Oftentimes, a large system, consisting of two or a few subsystems (degrees of freedom) 
interacting with their environment (thermal bath comprising a large number of degrees of freedom), can 
be adequately described as a composite system. Examples include electron transfer in solution 0, a large 
biological molecule, vibrational relaxation of molecules in solution, excitons in semiconductors coupled to 
acoustic or optical phonon modes. Quantum processes in condensed phases are usually studied by focusing 
on a small subset of degrees of freedom and considering the other degrees of freedom as a bath. 

Another interesting aspect is referred to the quantum to classical transition. The emergence of classicality 
is a typical consequence of having the quantum system in direct interaction with the external world. In fact, 
decoherence is the main ingredient in order to find classicality. The interaction between system and environ- 
ment induces a preferred basis which is stable against this interaction and becomes a classical basis in the 
Hilbert space of the coupled system. Preferred pointer states are resilient to the entangling interaction with 
the bath. This "einselection" (environment induced superselection) of the preferred set of resilient pointer 
states is the essence of the environment. It is accepted that a rapid loss of coherence caused by the coupling 
with the environment is at the root of the non-observation of quantum superpositions of macroscopically 
distinct quantum states 0. 

In this article, we analyze the decoherence induced in a composite quantum system, in which an observer 
can distinguish between two different subsystems, one of them coupled to an external environment. Our 
composite system is composed by a subsystem A coupled to a subsystem B which is also bilinearly coupled 
to an external environment £. The coupling to this external environment is only through subsystem B. 
Subsystem A remains isolated from £ but for the information delivered by B through a bilinearly coupling 
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between subsystems A and B. We will consider the thermal bath to be at high temperature and will work 
in the underdamped limit. 

In order to investigate this problem we mainly consider a simple model where subsystem A is represented 
by a harmonic oscillator and subsystem B is an upside-down one. The main motivation for studying this 
model is twofold. On the one hand, it is of interest to deepen and enlighten previous analysis of decoherence 
induced by chaotic environments. The upside-down oscillator has recently been used to model a chaotic 
environment which induces decoherence on the system 0- Even though it is an oversimplified model for a 
chaotic environment, it displays exponential sensitivity to perturbations, which is crucial in order to analyze 
chaotic evolutions. On the other hand, we want to emphasize that isolation from a chaotic environment is 
difficult, as has been noted in 0- Moreover, it is even harder to isolate a system from a chaotic environment 
than from the many harmonic oscillators of the quantum Brownian motion environment. In this context, 
we shall consider two different cases. Firstly, the case where the chaotic degree of freedom is part of the 
environment (i.e. an unstable system B) and is directly coupled to an external reservoir £ and to another 
subsystem A with different bare frequency. Secondly, the case where subsystem A is unstable and directly 
coupled to a harmonic oscillator (subsystem B) which is also coupled to an external bath £. These are 
the extension of previous works done in and 11, H for the first and second case, respectively. In both 
situations, we will estimate the decoherence time, which is the usual scale after which classicality emerges 
and is different for each case. We will show the dependency of these times with the parameters of the model. 

The analysis is completed by the inclusion of the other two different possibilities for the quantum composite 
system, i.e., a composite system constituted by a subsystem A coupled to subsystem B, both harmonic 
oscillators, and a composite system formed by subsystem A coupled to subsystem B, both inverted oscillators. 
As in the other two cases mentioned above, subsystem B is also coupled to an external reservoir £. All in 
all, we have four different composite systems to analyze. 

For every and each situation, we study the dynamics of the subsystem A. Not only did we study the 
influence of " its" environment (formed by subsystem B and £) at high temperature but also in the absence 
of the external reservoir £. In every case, we conclude that decoherence is faster in the case in which 
subsystem A is unstable. However, different conclusions can be arrived at when the subsystem A is a 
harmonic oscillator. All the cases studied in this paper have different decoherence time-scales associated, 
depending not only on the external temperature, but also on the type of subsystem one is considering in 
turn. Each case develops a different dynamics, being possible, sometimes, to find a quantum open system 
described using mixed quantum-classical dynamics [g (part of the composite system completely decohered 
while others didn't). 

In previous articles, a different composite system has been considered. For example, in Ref.^^l; subsystem 
A is taken to be a two-level system which is bilinearly coupled to a single harmonic oscillator B-subsystem; 
which is also coupled to an Ohmic (or supcr-Ohmic) set of infinite harmonic oscillators. Authors have shown 
that subsystem B losses coherence more rapidly than subsystem A, which maintains coherence for longer 
periods of time. This two-level composite system was also studied in in order to look at the exact 
solutions for the dynamics of the reduced density matrix of the composite system AB. Even though the 
composite system presented in this paper might appear to be similar to the one in [lol [Tlj , it has a different 
dynamics and therefore is of interest to investigate it separately. This article is an extended version of Ref. 

where we firstly presented our model and numerically evaluated the decoherence times in every case. 
Here, we are showing a complete analytical development of the influence functional method for the composite 
system (what is not included in From the influence action we evaluate the diffusive corrections to the 

master equation, and we complete the analysis about decoherence times, providing analytical and numerical 
estimations of the decoherence times, based on the inverted oscillators' dynamics in the phase space. 

The paper is organized as follows. In Section II we present our model and evaluate the influence functional 
for each of the considered cases, in order to compute, in Section HI, the diffusion coefficient of the master 
equation for the subsystem A. We evaluate diffusion analytically. Section IV is devoted to the analysis of the 
decoherence process. This is done by means of the decoherence factor. Section V contains our final remarks 
and in the Appendix we include details of the calculations. 
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II. THE MODEL: COMPOSITE QUANTUM SYSTEM IN AN EXTERNAL ENVIRONMENT 

A. General formulation 

We consider a ABE quantum system consisting of three coupled subsystems: subsystem A is coupled 
directly to subsystem B, while subsystem B is in direct contact with an external environment £. The total 
ABE classical action is, 

S[x,q,Q] = Sa[x] + SB[q] + SsiQ] + SAB[x,q] + Sssiq^Q]- (1) 

In the spirit of the quantum Brownian paradigm, the environment is taken to be a set of N independent 
harmonic oscillators with frequencies masses m„ and coordinates and conjugate momenta {Q,P) = 
{Qi, ...,Q]y, Pi, Pjy) so that the classical action is 

Se[Q]= fdsY^^iQl-ColQl). (2) 

Subsystem B consists of a single oscillator (upside-down or harmonic, depending on the case considered) 
with bare mass Mb, frequency fl and coordinate operator q, 

Sb{x)=^ j^ds^{e±^h^). (3) 
The interaction between subsystem B and the thermal environment is assumed to be bilinear, 

S^e^ I dsJ2cnq{s)Qn{s), (4) 
where c„ is the coupling constant to the nth oscillator. The environment is characterized by the spectral 

2 

density Is = tt^^ 2rral!r'^('^ ^ '^")- For simplicity, we assume an Ohmic environment, with the following 

spectral density l£{uj) — 2Af7o(I'e~ a^, where A is a physical cutoff, related to the maximum frequency 
present in the environment. 

Finally, we consider subsystem A consisting of a single oscillator (again, this oscillator can be an upside- 
down or harmonic one) with coordinate operator x whose classical action is 

Sa[x]^ f ds^^{x^±uj^x^). (5) 
Jo 2 

We suppose that subsystem A is bilinearly coupled to subsystem B by the interaction term 



^Xb = -A / dsxis)qis). (6) 







The dynamical properties of interest can be computed from the density matrix of the system at time t. 
The complete density matrix may be written in integral form in terms of the total propagator (we are setting 
the initial time to = 0) 

K{x,q,Q;t\xo,qo,Qo;0) = K{t\0) =< xqQ\cx]i{-iHt/h)\xoqnQo > 
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as 



p{x,q,Q,xi,qi,Qi) = J dxodxofdqodqQ/dQodQofK{t\0)p{0)k*{t\0). (7) 

We are primarily interested in the dynamics of the composite AB-system under the influence of the 
external environment £. In such case, the relevant quantity to analyze is the reduced density matrix pr, 
obtained by integrating out the environmental degrees of freedom. Such a reduction is especially correct if 
the characteristic time scale of the environment (which essentially is 1/A) is much shorter than those for 
the subsystem A and subsystem B. As is usually done, we assume a factorized initial condition between the 
composite system AB and the environment £, 

p{xo,xo/,qo,qo/,Qo,Qof;Q) = PAB{xQ,xo/,qQ,qo/;Q)p£{QQ,Qaf;0), (8) 

and the external environment initially in thermal equilibrium at temperature T. 

In this way we can write the integral form of the reduced density matrix at time t as 

Pr{x, xl, q, ql, t) = J dxodxoldqodqolJr{x, xl, q, ql; t\xo,xol, qo, qo'; 0)pab{xo,XoI, qo, Qo'; 0), (9) 
where the reduced time evolution operator Jr is 

Jr{x, xl, q, ql; t\xo,xol, qo,qo'; = ^ dQodQo'K{x, q, Q; t\xo, qo, Qo; 0)ps{Qo, Qo', 0) 

X K*{xl,ql,Ql;t\xal,qa',Qo';0)- (10) 

Given the initial conditions mentioned above, this expression for the reduced density matrix specifies a 
non-Markovian time evolution since the solution at time t depends on its past history. In the following we 
will use the influence functional method for deriving the master equation. This method provides us a way 
to obtain a functional representation of the evolution operator for the reduced density matrix. 



B. Influence functional method 



The formulation of the reduced density matrix in terms of an influence functional is widely discussed in 
the literature 

111 111 El- In the present paper, we will extend it to the composite AB-system, comprising 
two oscillators (harmonic or inverted). In the general setting, the evolution operator J for the full density 
matrix p is p(t) = J{t,0)p{0), where 

pxt rqt pQt 

J{xf,qf,Qf,Xi',qfl,Q'f\xo,qo,Qo,x'o,q'o,Q'o) = / Vx Vq PQe^'^^"^'?^^) 

•J Xn Jan Qo 

rQ't , , , 

Vq' / PQ'e-i^(^'«^'3). (11) 



rx{ 


Vx 


Mi 
/ 


1 Xq 




'go 




Vx' 


ri't 


lx'„ 







The path integrals are over all histories compatible with boundary conditions. As was mentioned in the 
Section above, our primary interest is in the effect of the external environment on our composite system AB. 
Therefore, we need the reduced density matrix for the AB system, defined as 

r+oo r + co 

Pr{x,x',q,q')= dQ dQ'p{x,q,Q\x',q',Q')8{Q-Q'), (12) 
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and the evolution in time is given by the reduced evohition operator Jr 

Pr{x,x' ,q,q';t) ^ / (ixodx^dgodgQ Jr{t\0)pAB{xo,x'Q,qo,q'o;0). (13) 



— OO J —OQ 



Assuming total separable initial conditions as mentioned above, the reduced propagator is 

r^i r^'i r1! pit 

Mxux'f,qi,q'f;t\xo,x'a,qo,q'o;0) ^ / Vx Vx' Vq p^'ei('SAW-SA(.)) 

J xo Jx'„ Jqo Jq'„ 

where F{x,x' ,q,q') is the Feynmann- Vernon influence functional JJ] (now for the composite system) given 

by 

dQi / dQn / dQ'f / VQ VQ' 

OO J — OO J — OO J Qq J Qq 

X ei(^^(Q)+^B.(9.Q)-S,(Q')-SB.(«',Q')) ps{Qo,Q'^) 

= eT;SA{x,x',q,q')^ (^5) 

where SA(x, a;', g, q') is the influence action for the AB composite system. Thus, we can define A{x, x' , q, q') 
as the coarse graining effective action: A{x,x' ,q,q') ~ Sa{x) — Sa{x') + SB{q) — Ssiq') + SAB{x,q) — 
Sab{x', q') + 5A{x, x',q, q'). 

It is important to note that in our model, the subsystem A is not directly coupled to the environment. 
Consequently, the influence functional is the well known influence functional 5 A for a bath of harmonic 
oscillators found in the literature (and only a function of q and q') 

SA{q,q')^-2f ds [ ds'Ag(s) r/(s - s') Sg(s') + » / dsf ds' Aq{s) v{s ~ s') Aq{s') (16) 
Jo Jo Jo Jo 

with 

Aq{s) = qis) q'is) ; Eg(,s) = l(g(.) + q'is)). (17) 

The kernels rj and v (dissipation and noise, respectively) are in general non-local and are defined as 77(5) = 
£l{s), with 



and 

K.) 



7(s) = / du) cosws, 

Jo ^ 

/ cttj if (tj) coth cosws, 

JO 2 



up to second order in the coupling constant with the external environment. In the high temperature limit, 
these kernels arc proportional to ~ 2AI^okBTS{s)/h and rj ^ AIjoS{s) 0,^^. As we are working in 
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the underdamped and high temperature hmit (which means fc^T ^ huj and 70 <^ ^2 but leaves ^okBT 
unrestricted), we can use the latter expressions. Therefore, if we evaluate Ea. H16|l . we have 

SA{q,qf) 2i -2AfB70 /* dsAq{s)tq{s) + i '^^^^^f^T f\^^^^^Yds. (18) 



" "'0 



Consequently, after integrating the bath of harmonic oscillators, we can write the influence functional 
F{q,q/), in the high temperature limit as 

F{q, qt) = e-"^ dsAa{s)tq{s)^^2M3^^ I^(Ag{.)fd.^ ^^g^ 
and therefore the redueed density matrix takes the form 

p,{x,x',q,q') ^ I dxodxo' dqodqo' dqfdqf/ Vq Pg/ 6^(^2(9) -Sb (9')) 



-00 J —00 "'—00 "'go "'go' 

Vxl e.i''SA{x)-SA{xl)) ^Ti\SAli{x,q)-SA'B(xl,ql)+SA[q,ql)) ^ 



J Xq J XqI 

At this stage we have all the information we need so as to estimate the effect of the thermal bath on the 
composite system AB. However, if we want to know how is the decohcrcnce process for the subsystem A, wc 
have to trace over all the degrees of freedom that belong to the new environment. If wc take a closed look 
at last expression, we notice that we can assume that our new problem is a subsystem A and a subsystem 
B which are coupled through an "effective interaction" Sl'^{x,q,x' ,q') defined by 

Sl'^{x,q,x\q') ^ 5AB(x,g) - 5AB(x^gO - 2MB7o£rfgAg(s)i]g(g) + ^ ^^'^^]_°^'^^ £dg(Ag(g))^ (21) 

C. The functional approach applied to subsystem A 

If we want to study the effect of subsystem B and the environment £ on subsystem A, we have to analyze 
the reduced density matrix but for this subsystem only. That is to say, we need 



00 /'OO 



Pr{x,xl)= j dq dq/p{x,q\x/,q/)5{q - q/), 
which is propagated in time by the reduced evolution operator J'r{x, xf) 

/OO poo 
dxo / dxofJj.{x,xf;t\xo,xof;0)pr{xo,xof;0). (22) 
-OO J —OO 

For simplicity we take that at t = the subsystem A and the new environment arc also uncorrclatcd, i.e., 
PAB{t = 0) = pA{t = 0) ® pBit = 0). Wc assume a Gaussian wave packet of the form e'^^'""®') ^/^'^ as 
the initial condition of subsystem B . This is a convenient choice in the sense that all of these states form a 
closed set under linear evolution 0, Then the evolution operator does not depend on the initial state of 
the system and can be written as in |l4| 

Jr{x{,Xf/;t\xo,xo';0) = Vx Da;/e^(^*(^)-^*(^'»JP(2;, a:/) = / Vx Vxf cxp I -A{x, xf) ^ , 

J Xq J Xq/ J XqI L J 

(23) 
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where we have defined T{x,xi) = e-k^-^i^--^') and A{x,x/) = Sa{x) — Sa{x/) +SA{x,x/) as the new influence 
functional and influence action, respectively. In order to evaluate 6A{x,xf), we must do the following 
integrations, 

/oo poo poc ^Qf pq^ 

dqo / dqo/ / dqi 'Dq Pg'e-fe(^«('?)-SB(g'))g^S'?a^,'?,^'.'?'), (24) 
- 00 J —00 J —00 J go J q'fi 

In order to perform the functional integrations, wc must solve the classical equation of motion for the 
subsystem B given by 



q{s)±n'q{s) = ^x{s). (25) 

In the latter expression, we have neglected the term related to the dissipation introduced by the external 
environment, because we assume an underdamped environment (small 70). 

At this stage, we must make clear that we want to analyze all possible combinations of subsystems A and 
B made up with a harmonic oscillator and one upside-down oscillator. Then, in some cases, subsystem B 
will be a harmonic oscillator (sign plus in Eq. (|25|) 'l and in others will be an upside-down oscillator (sign 
minus in Eq. (|25|l ). Wc will explicitly write the solution for one case, being possible to obtain the other 
solution by just replacing for ifl in the solution presented below. However, details are presented in the 
Appendix. 

Let's suppose subsystem B is an upside-down oscillator (case (a)) obeying q{s) — fl^q{s) = ■^x{s). In 
order to find the solution to this equation, we must find the solution to the homogeneous equation and to 
the particular one. After imposing initial and final conditions q{s = 0) — qo and q{s = t) ~ qi, respectively, 
we write the complete solution as 

sinh(r2(t - s)) sinh(17s) A sinh(f7s) /"*,., 
^^'^'^ ^ sinh(m) +'^^^hA(m)' A/B»sinh(m) X " ^))du 



A 



x(u)sinh(f](s - u))du. (26) 



Once the full expression for qc\{s) is known, we can go back to Ea. (|24|l and estimate it, obtaining as a 
result the effective infiuence action for subsystem A, 

5A{x,x/)=2 / dsi / ds2 y{si)i%si ~ S2)r{s2) + i / dsi / ds2 y{si)h'{si - S2)y{s2), (27) 

Jo JQ Jo Jo 

with y{s) = x{s) — xf{s) and r(s) = {x{s) + xf{s))/2. The quantities 77 and v may be defined as the new 
kernels of dissipation and noise, respectively, given by 

A^ 

^(si - S2) = ^^^^^ sinh(»(si - S2)), 

!>(si - S2) = ■^cosh(0(si - S2)). (28) 

In order to evaluate this new influence functional, we will use the saddle point method and, in this way, 
get rid of the functional integrals. Since the potentials in our model are harmonic, an exact evaluation of 
the path integral can be done. These integrals are dominated by the classical solution of the free equation 
of motion for subsystem A ^.t this stage, we assume that our subsystem A is a harmonic oscillator. 
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(being possible to obtain the solution for an upside-down oscillator by just replacing lo for icj), obeying 
x{s) + w^.t(s) = 0. 

If we ask for initial and final conditions of the form x{s = 0) = and x{s = t) = Xf . the classical solution 

is, 

, , sinfwft — s)) sin(ix's) 

{s) = xo \ ) " + X, — ^ . (29) 
sin(a;t) sin(a;tj 

Therefore, we can write the reduced evolution operator as in Ea. H23|) . and the reduced density matrix 
following Eq.lEHl. 

Finally, the reduced density matrix for the subsystem A takes the form, 

dxo dxo' Vx I Px' e^^(*)e-^(*VA(a;o,Xo;0), (30) 

with U and D, related to the unitary evolution and dccohcrencc process respectively, given by, 

/ , sin(w(t — s)) , , sin(a;s) „ /"* , [^^ , , ^ , ^ 

U={xo-xof) — \ \ ^, + [xi - xi>) . ) : - 270 / dsi ds2 2/(si)?7 Si - 52)^(52), (31) 
sm(wt) sin[ujt) ' ' 



Jo 



and 



D = f dsiAq,i{s)f + ^ / ds, I ds^y{s{)v{s, ~ S2)y{s2) (32) 



MT^ ./n ' ^ " 32ft 







From the last equation, we can see two contributions to the diffusion coefficient. The first one, proportional 
to the environmental temperature, comes from the coupling of subsystem B to the bath. The second term is 
the backreaction of subsystem B over A, through the weak A-coupling. Even though we are working in the 
high temperature limit, the underdamped bath (70 ^ ^2) produces both contributions could be of the same 
order of magnitude. Thus, both terms will be relevant in order to study decoherence effects on subsystem 
A. 



III. DIFFUSION COEFFICIENT IN THE MASTER EQUATION 

In this Section we will derive the diffusion coefficient in the master equation which will quantify the 
decoherence suffered by subsystem A for all four different cases. A commonly proposed way to analyze 
decoherence is by examining how the non-diagonal elements of the reduced density matrix evolve under the 
master equation. Following the same techniques used for the quantum Brownian motion to obtain the 
master equation we must compute the time derivative of the propagator JTI-, and eliminate the dependence 
on the initial conditions xq, xqI that enters through the classical solution Xc\{s). This can be easily done 
using the properties of the solution 



cos ^(t - s Af H i -ih—- 



AoJr(t, 0) 

where Aq — {xq — xqI), Af — (xt — Xfi) and Sf — (xf + Xfi) 



J,.(t,0) (33) 
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The master equation is commonly presented as 





" 92 


92 1 


[ 2A/a 







^M^n^ix^-x''] 



Pr{x,x';t) 



1 



+ -MASn^{t){x^ - x'^)p,{x,x']t) ~ ihr{t){x ~ x') 

- iMAV{t){x - x'fp,{x,x']t) - hT{t)f{t){x - x') 



' d__d_ 
dx dx 
d_ _d_' 

dx dx' 



Pr{x,x';t) 
p,{x,x';t), 



(34) 



where 'D(t){x—x')'^ is the diffusion term, which produces the decay of the off-diagonal elements. For simplicity 
we omitted the subindex f to indicate the final configuration Xf. Therefore, in order to find that coefhcient, 
we will take a closed look at those terms. The total diffusion coefficient is given by 



2jokBT .^2 j Agci(s) Aqci{s) + I ds v{t - s) Axci(s), 



a 
32h 



(35) 



where Aqc\{s) and i){t — s) are presented above in Eas. (|26() and H28() respectively. It is important to note 
that Agci(s) is the solution of the coupled system, and the new noise kernel is not the usual T-dependent 
noise kernel of the quantum Brownian motion. This term is solely coming from the interaction between the 
subsystems. 

In the rest of this Section we will present the exact results for the diffusion coefficients in all the four 
different situations considered. This can be summarized as: 



• Case (a): Harmonic Oscillator + Upside-Down Oscillator + £ 

This is the case we have developed so far of having a harmonic oscillator (subsystem A) coupled to 
an upside-down oscillator (subsystem B) by the interaction term presented in Eq. ©. The diffusion 
coefficient for this case is called 2?a- This is the generalization of the toy model considered in Ref.0 
where they did not consider the interaction of subsystem B (upside-down oscillator) with an external 
environment. It is easy to find results of just by setting 70 = in our results. We will plot this case 
in order to compare it with the other cases. Case (a) is the situation in which a Brownian particle (in 
a harmonic potential) suffers decoherence from an environment with one (or more) unstable degrees of 
freedom. 

• Case (b): Upside-Down Oscillator + Harmonic Oscillator + £ 

In this case we consider that subsystem A is an upside-down oscillator, obeying the classical equation of 
motion x{s) — lu'^x{s) ~ and subsystem B is a harmonic oscillator satisfying q{s) + Q'^q{s) = -p^x(s). 
It is straightforward to read the new solutions from the ones presented above by replacing lj ioj and 
ri — > ifi. The diffusion coefficient for this case is (see Appendix for details), and it represents the 
possibility of studying decoherence induced on an unstable system (toy model for a chaotic subsystem) 
by a completely harmonic environment 0, ITsL [T9| . We will see that this is the most decoherent system 
among all four cases studied in this paper. 

• Case (c): Harmonic Oscillator + Harmonic Oscillator + £ 

For completeness, we also consider the case of two harmonic oscillators coupled together by the in- 
teraction term Q and one of them (subsystem B) coupled to an external environment, very hot but 
underdamped. The procedure for deriving the diffusion coefficient is similar to what was done above 
but using the solution of a classical harmonic oscillator for subsystems A and B, as shown in the Ap- 
pendix. It is direct to find an analytical expression for this diffusion term (called 2?c) just by replacing 
— » iil in case (a). 
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• Case (d): Upside-Down Oscillator + Upside-Down Oscillator + E 

Finally, we consider two upside-down oscillators coupled together by the same interaction term than 
all other cases, and one of them coupled to E. This diffusion coefficient is Pd, and can be obtained 
from 2?b by making the change 17 — > ifi. We will see that this case is the most sensitive to external 
perturbations (both subsystems are unstable) when there is no external environment, thus decoherence 
is much effective than in the other cases. In particular, it is interesting to note that this case decoheres 
long before the others when there is no thermal environment (70 = 0). 

Once we have the analytical expression for the diffusion coefficients, we can study their behaviour for 
different relations between the parameters of the model. In order to illustrate some important cases, we 
present all four T) coefficients when both frequencies of the subsystems A and B are of the same order of 
magnitude (w « SX) and when w > fl, as shown in Fig. ^ Both cases are considered in the absence of 
external environment E (i.e. 70 = 0) and for low and high "foksT. Wc will restrict our results to times in 
which the no-dissipation approximation is still valid. 

In Fig. ^ we can appreciate the difference between the diffusion terms for composite systems when 
subsystem A is unstable and when it is not. In fact, the exponential behavior in the (b) and (d) cases is due 
to the fact that the final "decohered" subsystem A is unstable 0,0] (the upside-down osciUator). On the 
contrary, cases (a) and (c) present an oscillatory behaviour since subsystem A is a harmonic oscillator (and 
the solution Xci{s) for the classical equation of motion has oscillatory functions instead of hyperbolic function 
as in cases (b) and (d)). The difference between the exponential diffusion coefficients and the harmonic ones 
is particularly manifested for times cot > 1. For smaller times, all the coefficients are equivalent. But for 
longer times, cases (b) and (d) are easily dintinguishable from (a) and (c). The dynamics of the unstable 
upside-down oscillator is much in evidence after fit > 1, where cases (b) and (d) start to differenciate each 
other. 

The plots on top of Fig. are a generalization of the results obtained by Blume-Kohout and Zurek 
in which they considered a harmonic oscillator coupled to an inverted oscillator as an unique environment 
(there was no coupling between subsystem B and the external reservoir, i.e. 70 = in our model). Plots 
in the middle and at the bottom of the same figure are for small and high values of ^oksT since we have 
explicitly considered the presence of the external hot bath E, and thus, this presence modifies the diffusion 
terms adding a new contribution with respect to the one obtained in Ref. 0- 

In the case of 70 = 0, it is easy to note that (case (d)) grows faster (like an exponential function of 
time), while the other coefficients remain with a smaller amplitude. This is important in order to evaluate 
decoherence times, postponed till next Section. Instabilities inherent to the subsystem A exponentially 
enhance the diffusion originated due to the interaction with the environment. When the B-oscillator is 
also unstable, wc have more exponential sensitivity to perturbations than in any other case. However, it is 
important to note that it is a very toy model, in the sense that both oscillators are unbounded from below 
and therefore will develop un-physical oscillatory divergences j3- Therefore, wc conclude that in order to 
have a more precise idea of the physical consequences of having chaotic modes into the environment, wc can 
not neglect the interaction of these chaotic degrees of freedom with the rest of the world, modelled here as 
an infinite set of harmonic oscillators. 

The diffusion suffered by the subsystem A is the direct result of the interaction between A and B, and 
between the later and the external environment E when this last one is taken into account. The environment 
reacts to this interaction and the backreaction on the subsystem is by means of diffusion and dissipation. 
Let's take for example cases (b) and (d): unstable subsystem A is coupled to a harmonic oscillator (case (b)) 
or to an upside-down oscillator (case (d)). Subsystem A handles information to B via their coupling. In case 
(b), as B-oscillator is harmonic, the diffusion process is more effective. The B-oscillator has an oscillatory 
behaviour in time and is able of providing subsystem A with diffusion periodically. On the contrary, in 
case (d), the B-oscillator is unstable and unbounded from below. The stretching of its states is boundless 
(see next Section). Thus, part of the information is transferred to the reservoir, but, at large times, the 
intrinsic dynamics of the upside-down oscillator makes the provision of subsystem A with diffusion scarce 
and, consequently, less effective. A more quantitative explanation will be given in the following Section, 
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Case(b) ■ 
Case (c) 
Case (d) ■ 



Case (a) 
Case (b) ■ 
Case (c) 
Case (d) ■ 



Case(b) — 

Case(c) 

Case(d) ■■■■ 



Case (a) " 
Case (b) 

Case (c) 

Case (d) -■«■■■ 



Case (a) * 
Case (b) 

Case (c) 

Case(d) 



Case (a) 
Case (b) ■ 
Case (c) 
Case (d) ■ 



FIG. 1: We show the comparison between diffusion coefficients for tlie different cases considered. Plots on the top 
refers to the case 70 = 0, in which the composite system is isolated from the external bath. Plots in the middle and 
at the bottom represent the coefficients for 70 fcsT = 1 and 100, respectively. On the left column: we use uj — = 1, 
a = 0.01. On the right column: u — 5^1, a = 0.01. Cases (b) and (d) show a bigger growing rate at short times for 
both set of parameters. 



when we give an analytical estimation of the decoherence time for each case. 

Cases (a) and (c) are slightly different. Subsystem A is harmonic and the environment can have an 
unstable degree of freedom (case (a)) or not (case (c)). It is easy to see in Figure ^ that the backreaction 
of the full environment (B + £) on the subsystem A by means of the diffusion process is more effective for 
cases (b) and (d) (in this order) than for cases (c) and (a). This is due to the fact of having an inverted 
oscillator as the final subsystem A. In the middle of these figures, the "low temperature" limit is shown [we 
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are still working in the high temperature limit, but we are considering the underdamped case (70 is small 
with respect to any of the present frequencies), therefore the coefficient 'joksT could be a small number 
still preserving the hot bath assumption]. We can see that grows slightly faster than 2?d, and both of 
them are bigger than V^, and Vc- This can be understood by thinking in the dynamical properties of the 
composite system coupled to the external bath as we mentioned above. Oscillator B is not a good "diffusion 
handler" if it is unstable. At times longer than fit > 1, states in this oscillator are spread out too much. 
Diffusion must go from 5 to A through B (middle-environment subsystem), ft is not an effective process 
at short times (all the cases have a similar behaviour at short times). The dynamical behaviour of case (a) 
should be similar to the one occurring in case (b). However, the difference of having an inverted oscillator as 
the subsystem A or as the intermediate subsystem B is crucial. That is reflected in the exponential grow (or 
not) of the diffusion coefficient and in the decoherence times that we will present in the following Section. 

As wc stated before, in the (b) and (d) cases, the spreading of the initial state of the subsystem A is 
exponentially sensitive to fluctuations coming from the full environment (B + £), and it reacts quickly on 
the bath, losing information faster than in any other case. This is what happens in case (b), although the 
middle-environment is a harmonic oscillator. However, a question might arise: why I?b grows faster than 
fd, if the latter is "twice" unstable?. The key is: case (b) is the most decoherent because the unstable 
system losses its information in an uniform (non-ohmic) environment composed by a harmonic oscillator 
(B) plus an infinite set of harmonic oscillators in thermal equilibrium without internal unbounded regions 
(contrary to what happens in case (d)). In the low-T example, the major contribution to diffusion comes 
from the intrinsic dynamics of the composite system alone. At early times {tot < 1), unstable dynamics of 
the subsystem A dominates the temporal behaviour and, as they both have an inverted A oscillator, both 
cases (b) and (d) are similar (in fact, all the cases are similar at short temporal scales because at very short 
times both potentials are similar). However, when ujt > 1, and the presence of the external environment 
is still not so important, there is a noteworthy difference between I?b and I?d when the frequency of the 
subsystem A is similar to the one of B (see middle of Fig. QJon the left). When A has a bigger frequency, 
the dynamics is dominated by subsystem A and both diffusion coefficients are extremely similar in a longer 
temporal scale (middle Fig. fright). Pb and Pd arc indistinguishable ioi uj > il even at ^oksT ^ 1. 
However, as the temperature of the thermal bath increases, there is no distinction between cases (b) and (d) 
because the external environment dominates the diffusion coefficient. At high values of joksT we obtain a 
clear hierarchy between different composite systems (Fig. ^at the bottom). Again, it is easy to observe that 
the diffusion coefficients of those cases in which subsystem A is an upside-down oscillator reach bigger values 
than the others that have a harmonic oscillator as subsystem A. Therefore, we are able to conclude that the 
presence of instabilities into the composite system enhance decoherence. This effect is yet more important if 
the unstable subsystem A is coupled either to a solely chaotic degree of freedom (in the absence of external 
bath like our case (d)) or to an external environment (formed hy Ji + £) at high joksT, as in cases (b) and 
(d). In Fig. 121 we present the diffusion coefficient at a high value of ^^ksT for both systems which have a 
harmonic oscillator as subsystem A: cases (a) and (c). This figure should come in handy so as to compare 
these oscillatory coefficients with the hyperbolic-like other two. 



After integrating out all the degrees of freedom corresponding to the external hot environment Qn, and 
the coordinates q belonging to the subsystem B, we obtained the diffusive terms that induce decoherence 
on subsystem A. Therefore, we numerically integrated the diffusive terms in time, in order to plot the 
decoherence factor (see Appendix) 



Thus, Tit) is initially one (there is no interaction at f = between subsystems and environment), and it 



IV. DECOHERENCE IN A INTERACTING WITH B-£ 




(36) 
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Case (a) * 
Case (c) 



FIG. 2: Diffusion coefficients for cases (a) and (c) with uj — 5Q, a — 0.01 and ^oksT — 100. We can see the evolution 
of these coefficients at larger times than in Fig. 



decays to zero with time (this is the case of total decoherence). From the master equation for A-subsystem, 
it is easy to see that factor r{t) is at the root of the loss of quantum coherence. In order to illustrate the 
same cases, we present all four F coefficients when both frequencies of the subsystems A and B are of the 
same order of magnitude (uj w fl) and when w > f2, as shown in Fig. on the left and right columns, 
respectively. Both cases are considered in the absence of external environment £ (i.e. 70 = 0) and for low 
and high values oi^foksT. 

From the numerical results shown on top of Fig. O we can stress that in the absence of a hot bath, the 
decoherence time is smaller in case (d) than in (b), and both of them decohere long before cases (a) and 
(c). This is due to the fact that subsystem A. which is solely coupled to subsystem B, generates noise and 
dissipation at large scales. Thus, this noise and dissipation is bigger when the subsystem B is an upside-down 
oscillator (case (d)) than when it is a harmonic oscillator (case (b)). In this situation (70 = 0), case (d) is 
twofold exponential in time. In these figures, we can also observe what is going on for cases (a) and (c). In 
the former, the oscillatory dynamics of the A-oscillator and the hyperbolic stretching of the B-environment, 
proceed largely independently of one another. The B-environment induces only minor perturbations in the 
subsystem A and this subsystem does not disturb the environment. The stretching of the environment (due 
to being an inverted oscillator) along its unstable manifold is reflected in the system as diffusion. The same 
physical process occurs in case (b), with the sole and essential difference that the one stretching along an 
unstable direction is the subsystem A, while the environment is oscillating. As this stretching results in 
diffusion, the more stretching the system has, the more diffusion it feels. Case (d) is the best example 
in this "isolated" model, because both, A and B, stretch along a direction in the phase space, producing 
double exponential diffusion. This is the reason why it is the most decohcrent case. Case (c) is shown for 
completeness, but it is easily seen that decoherence occurs in a longer time scale (there is no stretching here). 

As soon as the interaction between B and the thermal environment is switched on, oscillator B dissipates 
not only on the bath but also on A. This is shown in the middle and at the bottom of Fig. |3| At very 
high environmental temperatures, there is no difference between cases (b) and (d); both of them decohere 
in the same temporal scale. The huge reservoir dominates the diffusion terms. But they still differ from the 
cases where there are harmonic oscillators as subsystems A (cases (a) and (c)). In Figure^ we present the 
behaviour of the T{t) coefficients for these cases for a longer time scale. We can observe that we need to 
wait longer times for decoherence to be effective in cases (a) or (c) with respect to (b) and (d) even in the 
highest temperature case. 
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FIG. 3: Decoherence factor for the same set than Fig. Isolated composite system decoheres first for the (d) case. 
For 7o 7^ 0, case (b) is more decoherent. 



A. Decoherence time predictions 



In this Section we will present an analytical estimation of the decoherence times based on the inverted 
oscillators' unstable dynamics of the phase space. 

When the final system A is an upside-down oscillator , an unstable point forms in the center of the phase 
space with associated stable and unstable directions. These are characterized by Lyapunov coefficients A 
with negative and positive real parts respectively. In order to have a quantitative expression for decoherence 
times we have to take into account that the dynamics now gives raise to the possibility of squeezing along 
the stable direction. The exponential stretching of the Gaussian packets in one of the directions due to the 
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FIG. 4: Cases (a) and (c) for 'yoknT — 100, ui = 50 = 5. Here we can see that, even the behaviour of r(t) is similar, 
case (c) decoheres first than case (a), in which B-system is an inverted oscillator. 



hyperbolic point is compensated by an exponential squeezing. The time dependence of the package width in 
the direction of the momenta is given by ap(t) = <7p{to) exp [At], where (Tp(to) is the corresponding width at 
the initial time. Lyapunov coefficient is given by the value corresponding to a linear potential A = 2cj^. 

Diffusion effec ts limit the squeezing of the Wigner function. The bound on the width of the packets is 
given by dc = \/2A7A H El (where i is b or d). There is another scale, tmax corresponding to the time 
in which decoherence starts to be effective, and after which squeezing becomes of the order of the limiting 
value. We use this to estimate the decoherence time scale. 

Evolution of the Gaussian packet will typically proceed in two different stages. During the first stage, 
evolution is dominated by the unitary part of the master equation and will be within an approximately 
preserved area. This lasts the time needed for the spreading of the initial state over a regular patch to be 
larger than the critical width. During this stage diffusion does not alter much the Wigner function, which 
is stretched or contracted. When the dimension of the patch becomes comparable with CTc, diffusion begins 
to dominate and the second stage of evolution begins. Further contraction will be halted at CTc but the 
stretching will proceed at the rate set by the positive Lyapunov exponent. As a result, the area (or the 
volume) in phase space increases. One can estimate the time corresponding to the transition from reversible 
to irreversible evolution as 



..41.^. (37) 

A CTc 

In our toy model, we can use this scale as the typical scale for decoherence, setting tjj ^ i_j^^p(^a^^ 
therefore we obtain 



is =imax+ yfil^^. (38) 

A (Jc 

For the same parameters used in Fig. 13 we are able to numerically estimate decoherence times as: 
toi, ^ 7.7 and ^ 6.4, for the first set of parameters on the left of Fig. ^{cu = fl = 1) where 70 = 0; 
toi, ^ 2.4 and to^ ~ 2.7; for joksT = 1, and tjj^^ ^ 1.6 and to^ ~ 1-7, in the hight T case joksT = 100. 
For the set on the right Fig. 13 (a; = 5fi = 5), we estimated: for 70 = 0; tof, ~ 3.0 and to^ ^ 2.7. We also 
got tuj^ J ~ 0.1, for "loksT = 1, and t^^^ ^ 0.6 in the case ^aksT = 100. All these results agree with the 
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decoherence times, defined by the times the decoherence factor T{t) goes to zero, tliat can be seen in the 
figures above. 

Using Ea. (|37|l . we can check that decoherence proceeds slower in case (d) than in case (b) for 70 ^ 0, 

to, - ^Dh = ^ In ^ = A- In — . (39) 
" A (T^ 2A Pb 

By simple inspection of diffusion coefficients in the figures, we can see that < Pb implying fjjj^ < to^- 

In the isolated from external environment case (70 = 0), we have 2?d > T^h resulting in t^^^ > tOd^ which 
agrees with our qualitative arguments and with what is shown in the plots on top of Fig. 13 

Decoherence times for cases (a) and (c) occur as for the usual harmonic systems. We can estimate them 
by using the result of the high temperature limit of the quantum Brownian motion paradigm, i.e. tu is the 
solution of: 1 sa J^^ 'D{s)ds (we have to take the typical distance L as 2a, proportional to the dispersion 
in position of our initial Gaussian packet). We present, in Fig. ^ T{t) for a longer time scale in order to 
establish the corresponding hierarchy in the environmental decoherent effectiveness. 



V. FINAL REMARKS 



In this article we analyzed the decoherence induced by an effective environment. The effective environment 
was considered to be formed by part of a composite system and an infinite set of harmonic oscillators. The 
composite system was considered to be any of the four possible combinations made up with a harmonic and 
an inverted oscillator. 

Since a set of harmonic oscillators is a stable system, small perturbations due to the state of the coupled 
system do not induce exploration of a large volume of the phase space for any oscillator. When one considers 
an inverted oscillator, it can explore its volume more efficiently when it is perturbed. 

When one works with a composite quantum open system in interaction with an stable environment, it is 
almost probed that the composite system (or part of it, for example its unstable degrees of freedom) will 
decohere before those systems where there are not inverted potentials (like our case (c)). Therefore, one 
could have a mixture of quantum-classical dynamics for the open composite system due to the fact that the 
different parts of the global system loss coherence at different rates. 

In our article we integrated out subsystem B, in order to study the effect of having (or not) unstable 
degrees of freedom into the full environment. Then we analyzed different situations and concluded that 
cases (b) and (d) are the most efficient (smaller decoherence times) at high temperatures, and (d) is the 
most diffusive case, when one turns off the thermal bath. There is a clear hierarchy between the different 
compositions of the composite systems. Those in which oscillator A is unstable (cases (b) and (d)) decohere 
before than those with a harmonic oscillator as the A-subsystem (cases (a) and (c)). At high temperatures 
of the external environment, it has been shown that cases (b) and (d) have the same decoherence time scale, 
while composite system (c) losses quantum coherence before case (a). 

As the system and environment interact, information about the initial state of the subsystem A is trans- 
ferred to the environment (and vice- verse). Case (b) is the most decoherent because the unstable system 
losses its information in an uniform (non-ohmic) environment composed by an harmonic oscillator (B) plus 
an infinite set of harmonic oscillators in thermal equilibrium without internal unstable (and unbounded) re- 
gions (contrary to what happens in case (d)). At intermediate temperatures, unstable A-subsystems decohere 
before stable A-subsystems because the intrinsic dynamics of those subsystems produce exponentially driven 
diffusive terms. However, at high temperatures the environment does not distinguish between subsystems 
and decoherence proceeds equally in every case. The effectiveness of diffusion depends on the intermediate- 
environment subsystem B. We have shown that harmonic oscillators keep information during a major period 
of time and therefore are able to react on the A-system more efficiently, than the case in which one has 
an inverted oscillator as "information delivery subsystem" . This is the main reason why case (b) is more 
decoherent (in general) than case (d) , and why case (c) is more effective losing quantum coherence than case 
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(a). It is important to stress that in order to have a more physical model for the effective environment which 
would contain an unstable degree of freedom one should consider a double well potential as B-oscillator. For 
example, for cases (a) or (d) . This non-linear potential has all the unstable properties of the upside-down 
oscillator at early times (giving a contribution to the diffusion terms which is exponential with time). Fur- 
thermore, this potential is bounded, which would be of much relevance in order to evaluate its global effect 
on subsystem A. In this situation, case (d) would be the most decohcrent case at any value of ^oksT [20l |. 

As was emphasized by authors of Ref . Q , the scale associated with the decoherence process when upside- 
down oscillators are taken into account as part of the coupled system, is logarithmically dependent on the 
coupling constant. This is easy to see from our analytical results, using the expression of the diffusion 
coefficient V, given in the Appendix fEa. (|A14|l ). into Ea. (|38|) . This implies that isolation from chaotic 
environments is "exponentially" difficult. It is even harder to isolate a system (or subsystem) from a chaotic 
environment than from the many harmonic oscillators of the quantum Brownian motion environment, where 
decoherence time is quadratic in the coupling constant. 
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APPENDIX A: DERIVATION OF THE DIFFUSION COEFFICIENT 



In this Section, we show the calculation of diffusion coefficient corresponding to case (a), in which we 
have a harmonic oscillator coupled to an upside-down oscillator, and it is in interaction with a set of infinite 
harmonic oscillators at temperature T. 

In order to perform the functional integrations of Eq. (|24() , we must solve the classical equation of motion 
for the subsystem B. If this system satisfies, g(s) — r2^g(s) = jg-x(s). The complete solution to this equation, 
after imposing initial and final conditions q{s = 0) = qo and q{s = t) — qi, is given in Ea. H26() . 

Once we have the classical solution for the upside-down oscillator coupled to a subsystem of coordinate 
x{s), we can write explicitly the influence action obtained after integrating out all the degrees of freedom of 
the environment £ fEa. (|18|l ). Neglecting the action of dissipation into the influence action (the dissipation 
coefHcient 70 is basically the square of the coupling constant between subsystem B and the environment £, 
[70 ~ c^] and we are working in the undcrdamped and high temperature limit (70 < fi < fc^T)), we can 
write it as 



t 

2 







dAiq,q/)^i2MBjoKT I ds(Ag(s))'. (Al) 

In last equation, Aq{s) is given by 



2 

sinh(m) J ' sinh(m) 



Aq.is) ^ (qo - qoff ( ^^^^^^ ) + 2(,„ - .0^) 5^«-^.9(., t) + ,(., tf , (A2) 



with 



g(s,t) = ^^(-[ duAx{u) sinh(»(s - u)) + / du Ax (u) smh(n(t - u)) ] (A3) 

MBil \ Jo smh{m) Jo 
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and Ax{u) = x{u) — xi(u). As the last term in expression (|A3|) does not depend on the initial conditions, 
it will be transparent to integrals in Eq. (|24|l . It is important to note that as subsystem B is an upside- 
down oscillator with frequency Vl, all the time-dependent functions in Eq. (|A3() arc hyperbolic trigonometric 
functions. 

After imposing pB{qo, 90'; 0) = ]\[g,^ii'}o-qo') )/2cr initial wave packet for subsystem B, the expression 

for the influence functional is 



f-l-00 n+00 p+00 i-qt nqf 

Hx.x') = c^p\UAiq,q')\= I dqJ dqo dqU Vq 2?<z'et*'^(«-*") (go, 9^,) 



I 

J cxD J ~OQ J ~oo J qo 'J q'Q 

X gi('S'B('?ol) + SAB(2:,gcl)-SB('3cl)-'SAB(a:',9d))^ ^^^^ 

where gci(s) is given in Ea. H26|l and 5A{q, qt) in Eq. (|A1() . Therefore, this integration can be done and yields 
(for a weakly coupled composite AB system) the result shown in Eas. (|27(l and (|28|l . 

Once the influence functional is known, it is straightforward to write down the reduced density matrix for 
subsystem A only 



p,{xf,Xf/;t) ^ dxo dxQi Vx P^'e^^-^-^^^^-^-^^^'^^e-f ('*.*)e^*-^(^'^'VA(a;o, xq/, 0), (A5) 
and the reduced evolution operator 

Jr{xi,X!r;t\xo,Xo/;Q)= Vx t^^,^USp^{^)-Sa{x,)) ^-g\s,t) ^l5A(x.x,) ^ (^g) 

In order to evaluate Ea. l)A6|l . we need the solution to the free equation of motion for subsystem A x[s) + 
uP'x{s) = 0. If we ask for initial and final conditions of the form x{s = 0) = and x[s = i) = Xf, the 
classical solution is, a;ci(s) — xt^^^^^^^^^y-^ + Xf , and the reduced evolution operator becomes 

J,{xf,Xii;t\xo,xoi;Q) = ei(SA{x^i)-SAix,^i))^-g\s,t)^j,SA{x,xO ^ ^i.u{t) ^-D(t) _ ^pj-^ 

with U and D, related to the unitary evolution and dccohcrencc process respectively, given by, 

I sin(w(t — s)) sin(ws) /"* f^^ , \~, \ / \ 

U={xo-xoi) . - ^. h {Xf - Xfi) . ~ 270 / dsi ds2y{si)ri{si - S2)r{s2), 

sm(wt) sm(wtj Jo ^0 

and 



_ 27oA^ 2 



ds(Ag(s))2 



32/1(7 



dsi 



ds2y{si)v{si - 52)2/(52)- 



Following the same techniques used for the quantum Brownian motion [l^ to obtain the master equation, 
we must compute the time derivative of the propagator Eq. IjA?!! . and eliminate the dependence on the 
initial conditions a;o, xqI that enters through the classical solution Xc\{s). This can be easily done using the 
properties of the solution 



AoJr(/,0) 



, , , , , sin((jj(/ — s)) ^ d 
cos{uj{t - s))Af -t- — ^-^ -T-^l^ 



Mt,o) 



(A8) 
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where Aq = [xq — xq/), Af = (xf — Xfi) and = (xf + Xf/). This identity allows us to remove the initial 
coordinates Xq, x'q by expressing them in terms of the final coordinates Xi, x[ and the derivatives dxf,dx'^, 
and obtain the master equation. 

The full equation is very complicated and, like in the case of the quantum Brownian motion, it depends on 
the system-environment coupling. In the present case, it also depends on the subsystems coupling constant 
A. As we are solely interested in dccohercnce, it is sufficient to calculate the correction to the usual unitary 
evolution coming from the noise kernel (imaginary part of the influence action). The result reads 

Pr{x[,Xf;t) - -i[i7i.cn,Pr] 
= -i[iJren,Pr] 

where ellipsis denote other terms that do not contribute to decoherence. 
This is equivalent to write 

Pr W -t[H,,n, Pr] ~ {xf " x't)^V{t)p„ (A9) 

with "D the diffusion coefficient fEa. (|A14|l ). Then, the effect of the diffusion coefficient on the decoherence 
process can be seen considering the following approximate solution to the master equation 



a_/27o 



kBT 



dt\ 



02 Jo 



ds g{sy 



dsds' Axc\{s)h'{s, s')Axc\{s') \pr 



Jo 



2g{s)g{s)ds + -^{xf ~ x'l) I ds cosh{n{t - s))Axci{s) ] pr 



where is the solution to the unitary part of the master equation ( i.e., without environment). The system 
will decohere when the non-diagonal elements of the reduced density matrix are much smaller than the 
diagonal ones. 

For the case of a harmonic oscillator coupled to an upside-down oscillator, which is coupled to an external 
environment, the diffusion coefficient fEa. (jA14p 'l for total the reduced density matrix p^{x[,Xf;t) is 



/ sinh(f2s) 



ds 



sinh(ris) 



n 



Vsinh(m) 
sin(a;(t — s)) + sin(wt) cosh(f2s) 



(cosh(f2i) cos(a;t) — 1) — cosh(51s) cos{ujt) + cos{uj{t — s)) 

sinh(r2s 



sinh(nt) 

cosh(rit)(l — cosli(fii) cos(u!t)) + sinhffJs) coshfwt)) + uj( — I s'miujt) coshfilf) 

/ V sinh(ilt) 



A^ 
32ft 



ds cos 



h{n{t ~ s))cos(w(t - s))|. 



(All) 



The procedure to evaluate the diffusion term for case (b), when one has an upside-down oscillator coupled 
to a harmonic one, and the last to the environment, is equivalent to what was done in the Subsection above 
but replacing ft ^ ifl and oj icj. This is because subsystem B satisfies q{s) + Q'^q{s) ~ j^x{s) and 

subsystem A satisfies x{s) — uj^x{s) = 0. 
Therefore, in this case, qci{s) is, 



9ci(s) 



sin(ri(t — s)) sin(r2s) A sin(Sls) 
qo . N 1- 9f 



sin(r2t) 



sm{nt) Mb^ sin(rjt) Jq 



.T(u)sin(ri(s — u))du 



x{u)sui{Q{s — u))du. 



(A12) 



20 



and Xci{s), 



Xcl[S) 



sinh(cj(t — s)) sinh(a;s) 
Xo r-TT — r^; \- Xf — 



sinh(wt) sinh(ti;t) 
If one follows all the procedure detailed above one obtains as the diffusion coefficient, 



(AI3) 



sin(r2s) 
sin(rit) 



{cos{ilt) cosh(a;t) — 1) — cos(f2s) cosh(ij;t) + cosh(w(t — s)) 



'.' ^„ cos(17t)(f — cos(f2t) cosh(ti>t)) + sin(rJs) cos(cji)^ + lo( — ^ | sinh(wi) cos(J7i) 



sm 



sm{nt)^ ^ ^ ' ^ " \ ) \ 'J \ sin(m) 

h([j(i — s)) + sinh(ci;t) 003(575)") + - ( dscoa{n{t — s))cos]\{Lo{t ~ s))\. 



(AI4) 



Diffusion terms I?c aud I?d can be obtained in a similar way, by replacing iQ in Pa and 2?b, 

respectively. 
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